Search for: All records

Histone modifications play a crucial role in regulating chromatin architecture and gene expression. Here we develop a multiscale model for incorporating methylation in our nucleosome-resolution physics-based chromatin model to investigate the mechanisms by which H3K9 and H3K27 trimethylation (H3K9me3 and H3K27me3) influence chromatin structure and gene regulation. We apply three types of energy terms for this purpose: short-range potentials are derived from all-atom molecular dynamics simulations of wildtype and methylated chromatosomes, which revealed subtle local changes; medium-range potentials are derived by incorporating contacts between HP1 and nucleosomes modified by H3K9me3, to incorporate experimental results of enhanced contacts for short chromatin fibers (12 nucleosomes); for long-range interactions we identify H3K9me3- and H3K27me3-associated contacts based on Hi-C maps with a machine learning approach. These combined multiscale effects can model methylation as a first approximation in our mesoscale chromatin model, and applications to gene systems offer new insights into the epigenetic regulation of genomes mediated by H3K9me3 and H3K27me3. 

Abstract The formation of condensed heterochromatin is critical for establishing cell-specific transcriptional programs. To reveal structural transitions underlying heterochromatin formation in maturing mouse rod photoreceptors, we apply cryo-electron microscopy (cryo-EM) tomography, AI-assisted denoising, and molecular modeling. We find that chromatin isolated from immature retina cells contains many closely apposed nucleosomes with extremely short or absent nucleosome linkers, which are inconsistent with the typical two-start zigzag chromatin folding. In mature retina cells, the fraction of short-linker nucleosomes is much lower, supporting stronger chromatin compaction. By cryo-EM-assisted nucleosome interaction capture, we observe that chromatin in immature retina is enriched with i ± 1 interactions, while chromatin in mature retina contains predominantly i ± 2 interactions typical of the two-start zigzag. By mesoscale modeling and computational simulation, we clarify that the unusually short linkers typical of immature retina are sufficient to inhibit the two-start zigzag and chromatin compaction by the interference of very short linkers with linker DNA stems. We propose that this short linker composition renders nucleosome arrays more open in immature retina and that, as the linker DNA length increases in mature retina, chromatin becomes globally condensed via tight zigzag folding. This mechanism may be broadly utilized to introduce higher chromatin folding entropy for epigenomic plasticity. 

Frameshifting is an essential mechanism employed by many viruses including coronaviruses to produce viral proteins from a compact RNA genome. It is facilitated by specific RNA folds in the frameshift element (FSE), which has emerged as an important therapeutic target. For SARS-CoV-2, a specific 3-stem pseudoknot has been identified to stimulate frameshifting. However, prior studies and our RNA-As-Graphs analysis coupled to chemical reactivity experiments revealed other folds, including a different pseudoknot. Although structural plasticity has been proposed to play a key role in frameshifting, paths between different FSE RNA folds have not been yet identified. Here, we capture atomic-level transition pathways between two key FSE pseudoknots by transition path sampling coupled to Markov State Modeling and our BOLAS free energy method. We reveal multiple transition paths within a heterogeneous, multihub conformational landscape. A shared folding mechanism involves RNA stem unpairing followed by a 5^′-chain end release. Significantly, this pseudoknot transition critically tunes the tension through the RNA spacer region and places the viral RNA in the narrow ribosomal channel. Our work further explains the role of the alternative pseudoknot in ribosomal pausing and clarifies why the experimentally captured pseudoknot is preferred for frameshifting. Our capturing of this large-scale transition of RNA secondary and tertiary structure highlights the complex pathways of biomolecules and the inherent multifarious aspects that viruses developed to ensure virulence and survival. This enhanced understanding of viral frameshifting also provides insights to target key transitions for therapeutic applications. Our methods are generally applicable to other large-scale biomolecular transitions. 

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases, including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition. 

The frameshifting RNA element (FSE) in coronaviruses (CoVs) regulates the programmed −1 ribosomal frameshift (−1 PRF) mechanism common to many viruses. The FSE is of particular interest as a promising drug candidate. Its associated pseudoknot or stem loop structure is thought to play a large role in frameshifting and thus viral protein production. To investigate the FSE structural evolution, we use our graph theory-based methods for representing RNA secondary structures in the RNA-As-Graphs (RAG) framework to calculate conformational landscapes of viral FSEs with increasing sequence lengths for representative 10 Alpha and 13 Beta-CoVs. By following length-dependent conformational changes, we show that FSE sequences encode many possible competing stems which in turn favor certain FSE topologies, including a variety of pseudoknots, stem loops, and junctions. We explain alternative competing stems and topological FSE changes by recurring patterns of mutations. At the same time, FSE topology robustness can be understood by shifted stems within different sequence contexts and base pair coevolution. We further propose that the topology changes reflected by length-dependent conformations contribute to tuning the frameshifting efficiency. Our work provides tools to analyze virus sequence/structure correlations, explains how sequence and FSE structure have evolved for CoVs, and provides insights into potential mutations for therapeutic applications against a broad spectrum of CoV FSEs by targeting key sequence/structural transitions. 

RNA motif classification is important for understanding structure/function connections and building phylogenetic relationships. Using our coarse-grained RNA-As-Graphs (RAG) representations, we identify recurrent dual graph motifs in experimentally solved RNA structures based on an improved search algorithm that finds and ranks independent RNA substructures. Our expanded list of 183 existing dual graph motifs reveals five common motifs found in transfer RNA, riboswitch, and ribosomal 5S RNA components. Moreover, we identify three motifs for available viral frameshifting RNA elements, suggesting a correlation between viral structural complexity and frameshifting efficiency. We further partition the RNA substructures into 1844 distinct submotifs, with pseudoknots and junctions retained intact. Common modules are internal loops and three-way junctions, and three submotifs are associated with riboswitches that bind nucleotides, ions, and signaling molecules. Together, our library of existing RNA motifs and submotifs adds to the growing universe of RNA modules, and provides a resource of structures and substructures for novel RNA design. 

Abstract The SARS-CoV-2 frameshifting element (FSE), a highly conserved mRNA region required for correct translation of viral polyproteins, defines an excellent therapeutic target against Covid-19. As discovered by our prior graph-theory analysis with SHAPE experiments, the FSE adopts a heterogeneous, length-dependent conformational landscape consisting of an assumed 3-stem H-type pseudoknot (graph motif 3_6), and two alternative motifs (3_3 and 3_5). Here, for the first time, we build and simulate, by microsecond molecular dynamics, 30 models for all three motifs plus motif-stabilizing mutants at different lengths. Our 3_6 pseudoknot systems, which agree with experimental structures, reveal interconvertible L and linear conformations likely related to ribosomal pausing and frameshifting. The 3_6 mutant inhibits this transformation and could hamper frameshifting. Our 3_3 systems exhibit length-dependent stem interactions that point to a potential transition pathway connecting the three motifs during ribosomal elongation. Together, our observations provide new insights into frameshifting mechanisms and anti-viral strategies.